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Insight into the large-signal operating characteristics of IMP A TT oscil- 
lators has been obtained from detailed numerical calculations of D. L. 
Scharfetter and H. K. Gummel. However, their calcidations are relatively 
expensive in computer time. Small-signal analyses of IMPATT diodes 
are much less expensive, but give little reliable infomiation about oscillator 
performance. 

This paper presents a model for large-signal analysis of IMPATT 
diodes which requires less than 1 percent of the computer time used by 
Scharfetter and GummeVs method, but still provides a realistic description 
of IMPATT oscillators in modes of operation which maintain carrier 
velocity saturation. We show graphical residts, based on numerical com- 
pidations, which provide information about phase relationships in 
IMPATT oscillators and improve understanding of the two-frequency 
mode of operation. 

I. INTRODUCTION 

Considerable insight into the large-signal operating characteristics 
of IMPATT oscillators has recently been obtained by D. L. Schar- 
fetter and H. K. Gummel through numerical calculations involving 
a complete modeling of the physical processes taking place in the 
diode. 1 They solve the one-dimensional partial differential equations 
for the generation, diffusion, and drift of holes and electrons, as well 
as Poisson's equation and the differential equations for the circuit 
in which the diode is embedded. Their analysis, which has led to 
improved understanding of IMPATT diodes and insight into new 
modes of oscillation, includes most of the important physical effects. 
Since their model is so detailed and accurate, their analysis is a 
large-scale computer project. Even after the computer programs are 
debugged, analysis of IMPATT diodes is relatively expensive. 

3S3 
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Another possible method is small-signal analysis, which is tractable 
numerically and more economical. 2-4 However, small-signal analysis 
gives only a limited amount of information about large-signal oscilla- 
tion. 

An intermediate model has been presented by W. J. Evans. 5 How- 
ever, his simplified large-signal model requires short transit-time in 
the drift region. He obtained a single first-order differential equation 
for the diode current, and was able to solve analytically for the 
(large-signal) diode current and voltage as a function of time. Unfor- 
tunately, the restriction to short transit time and his use of other 
approximations severely limit the applicability of his results. 

This paper presents a model for large-signal analysis of IMPATT 
diodes which incorporates most of the important physical processes 
in the diodes, but which is more tractable numerically than that of 
Scharfetter and Gummel. Because the new model requires far fewer 
computer hours than Scharfetter and Gummel's model, parameter 
studies are feasible; because the new model contains more physics 
than Evans' model or the small-signal models, the results are more 
reliable and more effects can be studied. 

II. LARGE-SIGNAL IMPATT MODEL 

We consider a one-dimensional diode, neglect diffusion, and assume 
that the electric field is large enough so that holes and electrons move 
at scatter-limited velocities. The latter assumption is essential for 
simplifying the calculations, but means that our model is not appli- 
cable to some important oscillator modes, such as analyzed by Schar- 
fetter, Bartelink, and Johnston. 6 

The other important simplification is that in Poisson's equation 
we neglect the space charge of the ac component of the particle 
current in the avalanche region. The detailed calculations of Schar- 
fetter and Gummel show that this is a good approximation for the 
classical Read mode of operation and other modes where the width 
of the avalanche region remains fairly constant during the oscillation 
cycle. 7 We include the space charge of the dc component of the par- 
ticle current in the avalanche region, and both ac and dc components 
of the particle current in the drift region. This says that the shape of 
the electric field profile in the avalanche region does not change in 
time, although the magnitude of the field does; that is, 

E„M, t) = 1(x) + E(t). (1) 
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This is equivalent to the commonly-made assumption that the particle 
current is constant in the avalanche region. 2 - 5 - " Given a bias current 
and the doping in the diode, we can calculate f(x) by solving a second- 
order differential equation by standard methods. 3 - 4 

We can treat diodes in which holes and electrons have unequal 
scatter-limitcd velocities and unequal ionization coefficients, although 
the model we have chosen for computation assumes the same values 
for electrons as for holes. The ionization coefficients may be any 
given function of the electric field, such as A e\p(—b/E). 

Besides the simplifications resulting from our diode model, we 
obtain important simplifications by solving directly for a periodic 
solution, rather than solving for a transient solution which eventually 
becomes periodic. We save computation time both because we elim- 
inate solving for the unneeded transient solution, and because we 
are able to use more advantageous numerical methods. 

Notice that the essential approximations are constant carrier veloci- 
ties and equation (1) ; many combinations of additional approxima- 
tions may be made, so that we are describing a whole hierarchy of 
IMPATT models. The aproximation of equation (1) can be over- 
come by a recursive method of solution, the first step of which is 
the subject of this paper. However, for most purposes equation (1) 
is sufficiently accurate. 

III. IDEALIZED DIODE 

In the rest of this paper, we consider a fairly simple model of an 
IMPATT diode; we are able to emphasize the essential features of 
the method without becoming lost in algebra. (We plan to treat the 
most general case in a subsequent article.) 

We consider a diode with a uniform electric field [/(.r) = 1] in the 
avalanche region (of length d a ) , and one drift region (of length d d ) . 
We assume equal hole and electron velocities and ionization coef- 
ficients, so that the continuity equations for holes and electrons in 
the avalanche region are 

dp(x, t)/dt = va[E(t)][p(x, t) + n(x, t)] - v dp(x, t)/dx 
dn(x, t)/dt = va[E(t)][p(x, t) + n(x, t)] + v dn(x, t)/dx. 

The importance of the approximations discussed in Section II is 
that they enable us to avoid dealing with the hole and electron den- 
sities as functions of both space and time. We proceed much as did 



386 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1969 

Read 7 and Evans, 5 adding the two continuity equations to obtain 
an equation involving the particle current density 

dl 9 (x, t)/dt = 2val v - qv 2 d(p - n)/dx. (2) 

The total current, particle plus displacement, is 

I At) = I v (x, t) + e BE(x, t)/dt. (3) 

I T is independent of x; since E in the avalanche region is independ- 
ent of x, I p {x, t) is also independent of x. Using this fact, we integrate 
equation (2) over the avalanche region. Neglecting reverse saturation 
current, we obtain 

dIJdt = ^{d a a[E(t)]-l\. 

If we use a (E) = A exp(—b/E), and define E c as the "critical" field 
(that field necessary to maintain the dc current under static condi- 
tions) then a(E ) = l/d a . We obtain 

dljdt = ^ [(Ad a y- E ' /Ei,) - 1]. (4) 

(Evans linearized the electric field, and used 



. 2v \E{t) "I 

dljdt =-■ -j- I v ml -jj£ - 1 J • 



d n 



We retain the more nonlinear version.) 

We treat the case in which there is an applied voltage* V(t) in ex- 
cess of that required to maintain the dc reverse bias current density 
I . The avalanche electric field is reduced by the space charge of the 
particle current in the drift region; as in Read, we have 

E(t) =E e + V(l)/(da + d d ) - \ J" 7,«')(l - •^ J ) dt' (5) 

where T d = d d /v is the time for particles to pass through the drift 
region after leaving the avalanche region. In equation (5) and the 
remainder of this paper we specialize to the case d a <K d d , as Read and 
Evans did. 

Equations (4) and (5), together with suitable initial conditions, 
may be solved for E(t) and I p (t). The total current, which is equal 



* We could equally well treat the case where there is an applied terminal 
current Ir(t). It is only slightly more difficult to embed the diode in a circuit; 
then we use Kirchoff's laws to write V(t) in terms of Ir(t) and the circuit 
parameters, generally as an integral over It. 
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to the terminal current, may be obtained from equation (3) : 

Ir(t) = j-j^j d dV(t)/dt + ^ f tTd U*) dt'. (6) 

The first term is the displacement current and the second is the par- 
ticle-induced terminal current, It p • 

Instead of solving the initial-value equations for the transient solu- 
tion, and extending the solution far enough in time for the limiting 
periodic solution to be obtained, it is preferable to solve directly for 
only the periodic solution. We assume V(t) to be periodic (not neces- 
sarily sinusoidal) with period t, and find solutions E(t) and I p (t) 
which also have period t. Now we need solve only for t in the inter- 
val (0, t) . It is more convenient to rewrite equation (4) as an integral 
equation; integrating equation (4) from to t, we obtain (T a = d a /v) 

I p {t) = 7,(0) exp ||- £ [{Ad a r E ' /E ^ - 1] d*'}. (7) 

We substitute equation (5) for E(t'), and have a single integral 
equation for I p (t), which we shall not write explicitly; since 
Ip(t + t) = I p (t), we solve only for ^ t ^ t. Because of its complexity, 
it will be solved numerically; for our method to be useful, we must be 
able to find numerical solutions quickly and reliably. We solve this non- 
linear integral equation by a standard iterative method, which converges 
quickly and stably. 

We divide the interval (0, t) into M intervals, let t n = nr/M, and 
solve equation (7) approximately for the M values I p (t n ), n = 1, 2, 
. . . , M. The integrals are evaluated by, for example, the trapezoidal 
rule, with integrands evaluated only at times t n . We obtain a system 
of M nonlinear equations for the values I p (t n ) , which are easily solved 
by, for example, the matrix analog of the Newton-Raphson method. 
More accurate integration rules than the trapezoidal rule may be 
used. They are more complicated to program, but save execution time 
because, for a given accuracy, fewer points need be used. 

After the I p (t n ) are found, then the Iritn) may be found from equa- 
tion (6). Then we Fourier-analyze I T and V; the ratio of their funda- 
mental components is the large-signal admittance at the frequency 
v = 1/t. The ac power delivered by the diode also is easily calculated. 

The small-signal admittance can be calculated analytically. It is 

Y( t ,\ = c [" 1 - cj'A)' 1 

W (d. + d d )T d U(«) + to/^rj 
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where 



*(«) =- 



exp (-icoT d ) - 1 + for, 

(-for,) 5 

and 



2 2a'(E c )vI 

0J o = 

e 

This admittance agrees with that of Gilden and Hines, when their 
"passive resistance of the inactive zone R 8 " is neglected. 2 

IV. RESULTS 

We analyze a germanium IMPATT oscillator similar to the one 
reported by Swan. 8 The results, which are presented graphically, are 
meant more as an illustration of typical results obtainable than as a 
comprehensive analysis of Swan's oscillator. In particular, the uniform 
avalanche region approximation is inadequate. However, even these 
preliminary results can clarify some aspects of IMPATT oscillator 
operation. 

The diode we analyzed has an avalanche width d a = 1.5 X 10 -4 cm, 
a drift width d d = 3.5 X 10~ 4 cm, and a carrier velocity v = 5 X 10 6 
cm per second. The ionization coefficient is &(E) = 5.9 X 10 +6 cm -1 
exp [-1.2 X 10 6 (V/cm)/E]. 

Figure 1 shows typical solutions for the ac components of V(t), I v {t), 
and I T „(t)- Two cycles are shown; we plot I Tp instead of I T , because the 
displacement current component of I T delivers no power. The frequency 
is 5 GHz, and the ac voltage amplitude is 10 volts. The avalanche 
region particle current is extremely nonsinusoidal, and becomes more so 
with increasing ac voltage amplitude. The fundamental component of 
I v (t) lags the ac voltage by about 49.5 degrees; this phase angle also 
increases with increasing ac voltage, approaching 90 degrees. The phase 
is closer to 90 degrees for higher frequencies; this may be attributed to 
the lessened effect of drift region space charge on the avalanche region 
field. (Read showed that the phase is 90 degrees when the space charge 
can be neglected, so that the electric field in the avalanche region is in 
phase with the terminal voltage.) The phase will be discussed further in 
connection with Fig. 3. 

Figure 2 is a complex-plane admittance plot; the dc component of 
I T is fixed at 500 A/cm 2 , and the ac component of V{t) is a pure sine 
wave, as it would be if the diode were embedded in a high-Q parallel 
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Fig. 1 — AC components of V, I P , 1t p at frequency 5 GHz, current density 
500 A/cm a , ac voltage amplitude 10 V. Two cycles shown. 



resistance-impedance-capacitance (RLC) circuit with only one res- 
onant frequency. The solid curves are curves of constant ac voltage 
amplitude; the dashed lines are curves of constant frequency. For 
frequencies about 6 GHz and above, the negative conductance de- 
creases with increasing ac voltage amplitude. For lower frequencies, 
the negative conductance first increases, then decreases with increas- 
ing ac voltage amplitude. For some low frequencies, for example 4 
GHz, the diode has positive small-signal conductance, but negative 
large-signal conductance for large enough ac voltage. This type of 
behavior agrees with that found at low frequencies by the detailed 
calculations of Scharfetter and Gummel. 9 

Insight into this type of behavior may be obtained from Fig. 3, which 
shows the amplitude and phase (with respect to the ac voltage) of the 
fundamental component of I p (t). Curves of constant ac voltage ampli- 
tude are solid and curves of constant frequency are dashed. To transform 
any point on this plot to a complex admittance plot, one first multiplies 
the complex number representing the point by 

[exp {-vaT d ) - l]/(-wjy 

to give the fundamental component of I Tp . (For small coT d , this is 
approximately the same as rotating the phase plot point of I v clockwise 
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200 




-20 

G in mhos per cm 2 

Fig. 2 — Complex-plane plot of diode admittance as a function of frequency 
and ac voltage amplitude. AC component of V(t) = A cos at. Current density 
500 A/cm 2 . 

by \u>T d .) If the resulting point is in the second or third quadrant, then 
the diode has negative conductance at that frequency and ac voltage 
amplitude. Then the amplitude is divided by the magnitude of the 
fundamental component of the ac voltage amplitude; finally the point 
is moved directly upwards by ew/(d a + d d ) to add in the displacement 
current component. 

At low frequencies, increased ac voltage amplitude improves the 
avalanche current phase considerably; at 4 GHz, where the small- 
signal conductance is positive, the improvement is sufficient so that 
the approximately %«>Ta rotation puts the phase plot point into the 
third quadrant, and negative conductance results for amplitudes 
above about 12 volts. 

Figure 4 shows output power in the fundamental, in watts per 
square centimeter. Maximum output power is obtained near the fre- 
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Fig. 3 — Magnitude and phase plot of fundamental component of avalanche 
particle current, l P (t), as a function of frequency and ac voltage amplitude. Cur- 
rent density 500 A/cm 2 . 



5000 
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Fig. 4 — Output power vs frequency for various ac voltage amplitudes. Cur- 
rent density 500 A/cm 2 . 
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quency with largest small-signal negative conductance. Efficiency 
depends upon the dc voltage across the diode, and therefore on the dc 
electric field in the drift region. The idealized model merely assumes 
this field large enough to maintain scatter-limited velocity. A better 
lower bound is obtained by calculating, at each point in the diode, 
the reduction in field resulting from the space charge of the current 
pulse drifting away from the avalanche region. Since It is independent 
of x, 

e dE(x, t)/dt = I p (t) + e dE aval (x, t)/dt - I v (t - x/v). 
Integrate from to t and subtract the dc component to obtain 

AE(x, t) = AJUitf) + f Uv(t') - hit' ~ x/v)] dt'. 

Jo 

Requiring the dc field at x to be as least as large as 

max [—AE(x, £)] 

gives the minimum field in the drift region allowable, and thus gives 
the minimum drift voltage and maximum efficiency. The maximum 
&E is also of interest, since in our model it is assumed that no ava- 
lanching takes place in the drift region. If E a is the minimum electric 
field for which appreciable avalanche ionization can take place, then 
we must require that 

max [AE(x, t)] + max [-AE(x, t)] < E a 

for all x. This limits the allowable ac voltage amplitude. Numerical 
results indicate that, for a given current density, the allowable ac 
voltage amplitude depends only weakly on the frequency of oscillation 
and the shape of the voltage waveform V(t). For our sample diode, if 
E a « 10 5 volt/cm, 20 volts ac amplitude is, approximately, the limiting 
voltage for the simple diode model used. A more general model, with a 
nonuniform avalanche region, could partially include avalanche- 
broadening, and thus be used to analyze diode modes of operation with 
larger ac voltage amplitudes. 

Swan, Lee, and Standley recently reported marked improvement 
in the performance of an IMPATT oscillator embedded in a micro- 
wave circuit resonant at two frequencies, the fundamental and its 
second harmonic. 8 - 10 As an aid in understanding their results, in the 
rest of this section we present calculations in which V(t) has the 
form 7 cos «£ + %V sin 2oit. (This particular waveform is not the 
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optimum form, but is illustrative.) "Best" values of the magnitude 
and phase of the second harmonic depend on the frequency and mag- 
nitude of the fundamental, the bias current, the particular diode 
used, and the standards of performance to be applied to the diode. 

Figures 5 to 8 correspond to Figs. 1 to 4; V is now chosen so that 
half the peak-to-peak value of V(t) is 5, 10, 15, or 20 volts. In com- 
paring the oscillations resulting from a two-frequency V(t) with a 
one-frequency V(t), we have chosen to compare voltage waveforms 
with equal ac peak-to-peak amplitudes, not with equal amplitudes 
in the fundamental, because of the ac amplitude limitation discussed 
earlier. 

Figure 5 shows solutions for the ac components of V(i), I v (t), and 
I T Jt)', the frequency of the fundamental is 5 GHz, and the maximum ac 
voltage is 10 volts. The fundamental component of I p (t) lags the funda- 
mental component of the ac voltage by about 60 degrees instead of 
49.5 degrees. 

Figure 6 is a complex-plane admittance plot; for clarity, the lines 
of constant ac voltage are omitted, and dots are placed at 0, 5, 10, 15, 
and 20 volts. For nonsinusoidal ac voltages, the admittance at the 
fundamental can move to the left of the small-signal admittance 
line, resulting in a better negative conductance-to-susceptance ratio. 
This is especially apparent at 4 GHz; comparison of the 4 GHz, 20 




TIME 



Fig. 5 — AC components of V, I p , I Tp at frequency 5 GHz, current density 500 
A/cm 2 , ac voltage amplitude 10 V. Two cycles shown. 
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Fig. 6 — Complex-plane plot of diode admittance at fundamental frequency 
as a function of frequency and maximum ac voltage. AC component of V(t) = 
Vocos at + iVoSin 2ut. Current density 500 A/cm a . 



V points in Figs. 2 and 6 indicates that the addition of a second 
harmonic component has improved the diode negative conductance 
by a factor of three, and the negative conductance-to-susceptance 
ratio by a factor of four. 

Figure 7, corresponding to Fig. 3, shows in more detail the phase 
improvement obtained with addition of the harmonic component to 
the ac voltage. 

Figure 8 shows output power in the fundamental. At higher fre- 
quencies, maximum output power is less than that shown in Fig. 4. 
At 5 GHz, maximum output power is about the same as in Fig. 4, 
but occurs at an improved negative conductance-to-susceptance ratio. 
This ratio is important when the diode's parasitic resistance R p , which 
we have neglected, is included, since the power is proportional to 
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Fig. 7 — Magnitude and phase plot of fundamental component of avalanche 
particle current, 7,(4), as a function of frequency and ac voltage amplitude. 
Current density 500 A/cm 2 . 

— Q/(G 2 + B-) — R p . At 4 GHz, maximum output power is sub- 
stantially better than that shown in Fig. 4. The nonsinusoidal form 
of V(t) improves the phase of l T {t) sufficiently so that, even though 
V(t) contains less of the fundamental, the output power in the fun- 
damental remains about the same. 
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Fig. 8 — Output power in fundamental vs frequency for various maximum 
ac voltages. Current density 500 A/cm 2 . 
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V. CONCLUSIONS 

We have presented a numerically tractable method for analyzing 
large-signal IMPATT oscillators. It requires less than 1/100 of the 
computer time used by Scharfetter and Gummel's method, but still 
provides a realistic description of IMPATT oscillators in modes of 
operation which maintain carrier velocity saturation. (One solution 
takes about one second of computation time on the GE 645 computer, 
if we use 15 time steps per period.) Numerical calculations based on 
a simplified version of the model have provided information about 
the phase relationships in an IMPATT oscillator and improved un- 
derstanding of the two-frequency mode of operation. 
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